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Abstract. The effect of density dependent symmetry energy is studied on elliptical flow is studied us- 
ing isospin-dependent quantum molecular dynamics model (IQMD). We have used the reduced isospin- 
dependent cross-section with hard(H) equation of state to study the sensitivity of elliptical flow towards 
symmetry energy in the energy range of 50 - 1000 MeV/nucleon. The elliptical flow becomes zero at a par- 
ticular energy termed as transition energy. A systematic effort has been made to pin down the transition 
energy for the density dependent symmetry energy. 
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1 Introduction 

Information about the nuclear matter under the extreme 
conditions of temperature and density and the role of sym- 
metry energy under these conditions is still a topic of 
crucial importance in the present day nuclear physics re- 
search. Quite good progress has been made with the help 
of heavy-ion collisions to study the nature of equation of 
state [Il[2]- The form and strength of symmetry energy, 
and the study of different observables at sub-nuclear den- 
sities with respect to the symmetry energy is still a chal- 
lenging task in the intermediate energy heavy ion colli- 
sions(HIC's). Among different observables, collective flow 
and its various forms enjoys a special status, due to its sen- 
sitivity towards the model ingredients that define a equa- 
tion of state. 

Recently, many efforts have been made to study the 
collective flow both theoretically and experimentally in 
heavy-ion collisions [glll[51|gl[71[gl|gl[n] . The collective flow 
is a motion characterized by the space-momentum corre- 
lations of dynamic origin. The directed flow and elliptical 
flow are the two parameters available in the heavy-ion col- 
lisions to study the collective flow. The directed flow seems 
to be constrained only along the reaction plane which is 
due to the bounce-off of the compressed matter. However, 
the elliptical flow is more suited to study the collective 
flow which is a squeeze out of the spectator matter out of 
the reaction plane [51 ll0l[TT| . The highly stopped and com- 
pressed nuclear matter around the mid-rapidity region is 
seen directly in the squeeze out. Generally, the elliptical 
flow is the difference between the major and minor axes. 
It describes the eccentricity of an ellipse like distribution. 



The elliptical flow is defined by the second-order Fourier 
coefficient from the azimuthal distribution of detected par- 
ticles at mid-rapidity. 



Mathematically, 



dN 



—— = Poil + '2vicos(j) + 2v2Cos2(l)) . 



(1) 
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Here, 4> is the azimuthal angle between the transverse 
momentum of the particle and reaction plane. The param- 
eters {Cos2(j)) of elliptical flow depends upon the complex 
interplay between the expansion, rotation and shadowing 
of the spectators, apart from the incident energy. The par- 
ticular energy at which elliptical flow vanishes is termed 
as the transition energy. The elliptical flow arises due to 
the orthogonal asymmetry in the configuration space and 
rescattering. 

One needs to understand that the high density symme- 
try energy term can be probed from the isospin effects in 
HIC's. The elliptical flow, which provides us an opportu- 
nity to study the pressure that is generated very early dur- 
ing the reaction [10], is highly sensitive towards the sym- 
metry energv[9. .12tll3lll4] . The studies based on the Trans- 
port model calculations concluded, that the squeezed out 
of nucleons with high transverse momenta perpendicular 
to the reaction plane, i.e. elliptical flow, has probably the 
highest sensitivity towards the symmetry energy [Tl]. 

Recently, an attempt was made by S. Kumar el. al.[S] 
to correlate the elliptical flow with isospin content of a 
reaction. Various studies [91IT3] concluded the elliptical 
flow as a powerful probe for the symmetry energy. Unfor- 
tunately, all the calculations were silent about the den- 
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sity dependence of the symmetry energy. The symme- 
try energy which is the difference of the energy per nu- 
cleon between pure neutron matter and symmetric neu- 
tron matter is taken as 32 MeV corresponding to the nor- 
mal nuclear matter density i.e. p = 0.16 fm~'^ P^f^U] . This 
understanding does not remain valid as one goes away 
from the normal nuclear matter density and symmetric 
nuclear matter [14J[2D]. The knowledge of the density de- 
pendence of the symmetry energy and its equation of state 
for isospin asymmetric nuclear matter is of crucial impor- 
tance, to study the structure of the systems as diverse 
as the neutron rich nuclei and the neutron stars. A large 
number of studies have been performed on the density de- 
pendence of the symmetry energy in the recent past [HI 
[T51I16) . The equation below gives us the most extensively 
used theoretical parametrization of how the symmetry en- 
ergy varies against p P^[T31ITB] . 



enters into the calculations via symmetry potential, cross- 
sections and Coulomb interactions [21]. The details about 
the elastic and inelastic cross-sections for proton-proton 
and neutron- neutron collisions can be found in Ref.|21j. 

In IQMD model, the nucleons of target and projectile 
interact via two and three-body Skyrme forces, Yukawa 
potential and Coulomb interactions. In addition to the 
use of explicit charge states of all baryons and mesons, a 
symmetry potential between protons and neutrons corre- 
sponding to the Bethe- Weizsacker mass formula has been 
included. 

The hadrons propagate using classical Hamilton equations 
of motion: 



drj 
dt 



d{ H 
dpi 



dpi 
dt 



d{ H 

dn 



(3) 



with 



E{p) = E{p,){p/p,)\ 



(2) 



The larger(smaller) value of the constant 7 corresponds 
to stiff(soft) density dependence of the symmetry energv[T^ 
I15l[l6, . Experimentally, symmetry energy is not a directly 
measurable quantity and has to be extracted from the ob- 
servables which can shed light on symmetry energy. Some 
authors from the FOPI data concluded that the symmetry 
energy at supra-densities is very soft [T7]. However, oth- 
ers concluded an opposite findings [18]. Recently, Shetty 
et.al.[Tn], and Tsang et. al.[Tn], concluded the density de- 
pendence of the symmetry energy with 7 = 0.4 - 1.05. 

In a recent communication, it was advocated that at 
high baryon densities, large symmetry energy repulsion 
affects the transverse in-plane flow [5D] of a neutron rich 
matter in heavy-ion collisions. A large sensitivity for the 
transverse in-plane flow was detected by assuming the 
symmetry potential as a function of density. However, it 
could be of the interest to carry out a meaningful investi- 
gation for elliptical flow, by including the density depen- 
dence of the symmetry energy. No such study is reported 
in the literature so far. Our present aim is to pin down the 
density dependent symmetry energy via transition energy. 

The present study is carried out within the frame- 
work of the isospin-dependent quantum molecular dynam- 
ics(IQMD) model [5T]. Section II describes the model in 
brief. Section HI discusses the results, and Section IV sum- 
marizes the results. 



2 ISOSPIN-dependent QUANTUM 
MOLECULAR DYNAMICS (IQMD) MODEL 

The IQMD model [5T], which is an improved version of 
the QMD model [^[^ developed by J. Aichelin and 
coworkers and then improved by Puri and coworkers has 
been applied to explain various phenomenon such as col- 
lective flow, disappearance of flow, fragmentation & el- 
liptical flow successfully. The isospin degree of freedom 



{H)^{T)+{V) 

„2 



^ 2m,- 



EE 

I j>i 



Mr,p,t)V'^{r',r) 



xfj{r',p', t)drdr' dpdp' . 



(4) 



The baryon-baryon potential V''^ , in the above relation, 
reads as: 

V^^ (r — r) = Vg-^y^j^^ + Vy^kawa + ^Coul + Kym + "^mdi 

= (hd{r' ~r)+t2d{r' ~r)p^-^ 



+ t: 



exp{\r' — r\/p) 
' i\r'-r\/p) 



\r' — r\ 



+t,-T^TiS{r/ -rj) 
+t7ln^[ts{pi' -pf + l]S{rl-r) 



(5) 



■th 



Here Zi and Zj denote the charges of i and j 
baryon, and TJ, Tg are their respective T3 components (i.e. 
1/2 for protons and -1/2 for neutrons). The parameters p 

and ti, , are adjusted to the real part of the nucleonic 

optical potential. For the density dependence of nucleon 
optical potential, standard Skyrme-type parameterization 
is employed. The potential part resulting from the convo- 
lution of the distribution function with the Skyrme inter- 
actions Vskyrme fcads as : 

^Skyrine — 



Pint 
PO 



P f Pint 
\ PO 



(6) 



The two of the three parameters of equation of state are 
determined by demanding that at normal nuclear matter 
density, the binding energy should be equal to 16 MeV. 



9p^ 



dp^ \A 



E 



(7) 



K is the compressibility factor. The different values of com- 
pressibility give rise to soft(S) and hard(H) equations of 
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state(EOS). It is worth mentioning that as shown by Puri 
and coworkers, Skyrme forces are very successful in the 
analysis of low energy phenomena such as fusion, fission 
and cluster-radioactivity, where nuclear potential plays an 
important role. |24j . The momentum dependent interac- 
tions can be incorporated by parametrizing the momen- 
tum dependence of real part of optical potential[25j. 

Many studies concluded the matter to be soft, whereas 
many more believe the matter to be harder in nature [261 
[27] . As noted ^28;, elliptical flow is unaffected by the choice 
of equation of state. For the present analysis, a hard (H) 
equation of state has been employed along with reduced 
isospin-dependent cross-section(0.9 of <jnn)- 



3 Results and Discussion 



We here perform a complete systematic study for the reac- 
tion of Aug"^^ + Augj^, using different parametrizations of 
the density dependence of symmetry energy. In addition to 
that, the reactions oi SnU'^' + Snll'^ and Xe^^^+Xe\\^ are 
also simulated to acknowledge the effect of isospin content 
on the elliptical flow. The phase apace generated by the 
IQMD model has been subjected to clusterization using 
minimum spanning tree method[MST] |24) and analysis 
packages [23]. The MST method binds two nucleons in a 
fragment if their distance is less than 4 fm. However, in 
recent times many effective and proficient, algorithms are 
also available in the literature. The quantum molecular 
dynamics(QMD) generated phase space analysed through 
one such algorithm simulated annealing clusterization al- 
gorithm (SAC A) [29] tends to solve the long standing prob- 
lem of QMD method. The reactions are followed until the 
elliptical flow saturates i.e. at 200 fm/c. 



To study the effect of symmetry energy. We display in 
fig. 1, the reactions of ATe^f + Ae^f and Snlf + Snlf 
under the same conditions for LCP's. Interestingly, the el- 
liptical flow seems to be more sensitive for the N/Z ratio 
at incident energy of 100 MeV/nucleon. The trends ob- 
served through our calculations shows that irrespective of 
the same mass, the weaker squeeze-out flow is observed in 
case of the neutron rich system. These findings are also 
in accordance with the findings in ref. [9,J50 . At the inci- 
dent energy of 50 MeV/nucleon, the isospin does not play 
a larger role. In such case, binary collisions are rare and 
mean field dominates the dynamics of the reaction. On the 
other hand, around 100 MeV/nucleon, the mean field and 
binary collisions both contribute towards the dynamics of 
the reaction. In such a case larger effect of isospin can 
be observed. Therefore, it could be further interesting to 
carry out the study of elliptical fiow for the various forms 
of the density dependence of the symmetry energy. 
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Fig. 1. Transverse momentum dependence of elliptical flow, 
summed over entire rapidity distribution, for LCP's at 50(top) 
and 100 MeV/nucleon(bottom), respectively. Both reactions 
have same mass number but different N/Z ratio. 



0.4 



0.2 



0.0 




A 

> 0.2 



0.0 
0.3 

0.2 

0.1 

0.0, 




197 . , 197 , 

Au,„ + Au,„ 1 .0 



79 79 

/\e = 100 AMeV 




0.5 



0.0 
-0.4 




0.2 

0.0 
0.2 

0.1 

0.0 



100 200 300 400 100 200 300 400 500 

< Pj>(MeV/c) 

Fig. 2. Transverse momentum dependence of the elliptical 
flow, summed over the entire rapidity distribution, at b = 5 
fm for the system Aii^-^ -\- Aiigy^ for different forms of the den- 
sity dependence of the symmetry energy, at 50 (left) and 100 
(right) MeV/nucleon. The top, middle, and bottom panels rep- 
resent the free nucleons (FNs), light charged particles (LCPs), 
and intermediate mass fragments(IMF's) respectively. 



As the symmetry energy acts directly on the the squeeze 
out nucleons, that are mostly from the high density region, 
formed during an early stage of the reaction and are not 
affected by the spectator nucleons [14]. For neutron rich 
Au^-f' large role of symmetry energy can 
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Fig. 3. Transverse momentum dependence of the elliptical 
flow, for mid rapidity (| y |= "'■ < 0.1), at b = 5 fm for 
the system Aul^^ + Au^'^ for different forms of the density de- 
pendence of the symmetry energy, at 100 (right) MeV/nucleon. 
The top, middle, and bottom panels represent the free nucleons 
(FNs), light charged particles (LCPs), and intermediate mass 
fragments(IMF's) respectively. 



be expected with reference to the elliptical flow. The study 
of elliptical flow by constraining the density dependence 
of the symmetry energy can give us better understanding 
of nucleon-nucleon interaction and its extrapolation to the 
behaviour of symmetry energy in neutron-rich nuclei. In 
fig. 2, we display the, final state elliptical flow, at different 
incident energies, for the free particles (upper panel), light 
charged particles (LCP's) [2 < A < 4] (middle) and inter- 
mediate mass fragments [IMF's] {5 < A < Afoj /6] (lower 
panel) as a function of transverse momentum (Pt ) for the 
various forms of the density dependence of the symmetry 
energy i.e. 7 — 0.66, 1.33 and 2, respectively. The pressure 
created during a collision can be revealed through the el- 
liptical flow, with respect to transverse momenta. 




1.5 1.0 1.5 2.0 100 120 140 160 
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Fig. 4. Same as in fig. 3, but for the transition energy depen- 
dence of gamma(7). 



flow on a large scale. 

The isospin effect basically generates from the mid- 
rapidity region or in other words from the participant 
zone. To study the effect of symmetry energy, we parametrized 
the {V2) as a function of (Pt) at mid-rapidity (|2/| = — 
0.1) in fig. 3. We display the FN's(top panel), LCP's('nhddle 
panel) and IMF's(bottom panel) for the various forms 
of the density dependence of symmetry. We see that the 
squeeze out is affected more in case of LCP's as compared 
to other fragments. Although, the heavier fragments have 
weaker senstivity towards the symmetry energy. 

The different forms of the symmetry energy tends to 
affect the elliptical flow considerably. The most interesting 
aspect in this case can be the variation of the transition 
energy (ETrans) with different forms of the density de- 
pendent symmetry energy. The incident energy at which 
elliptical flow vanishes is termed as the transition energy. 
Although, the variation in the squeeze out of nuclear mat- 
ter due to the various forms of the symmetry energy can 
be understood with the transition energy. 



This elliptical flow is integrated over the entire rapid- 
ity range. Gaussian-type, behaviour is observed for all the 
forms of the density dependence of the symmetry energy 
i.e. for the different values of gamma. 

The effect of symmetry energy is clearly visible in the 
figure. The elliptical flow seems to be sensitive towards 
the various forms of the density dependence of the symme- 
try energy. The trends observed through our simulations 
shows the weaker squeeze out flow for the larger values of 
gamma. The very stiff form of symmetry energy i.e. 7 = 2, 
corresponds to the maximum positive value of elliptical 
flow. Indeed, the very stiff form of symmetry energy i.e. 
7 — 1.33 and 2 does not seems to affect the squeeze out 



In fig. 4, we parametrized the transition energy with the 
stiffness factor 7. The transition energy tends to increase 
with the increase in value of gamma. The light charged 
particles however tends to behave in a different manner 
as compared to the intermediate mass fragments. As, the 
FN's and LCP's originate from the participant zone. The 
larger sensitivity of free nucleons and light charged par- 
ticles tends to increase the transition energy upto a cer- 
tain level. The (ETrans) tends to decrease for very stiff 
form of the symmetry energy. Although, the IMF's are 
being generated solely from the spectator zone. The tran- 
sition energy for the IMF's increases with 7. As the IMF's 
being heavy particles move slowly as compared to FN's 
and LCP's. The symmetry energy tends to play a role 
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Fig. 5. Variation of the elliptical flow, summed over the entire 
transverse momentum, with beam energy at \y\ — < 0.1 

for the reaction Au^^ + Awg™ . Here a comparison is shown 
with the constant form of symmetry energy and the density de- 
pendent symmetry energy for gamma = 0.66, with the exper- 
imental flndings of the INDRA, FOPI, and PLASTIC BALL 
Collaborations gllSlle] 



at later stages in case of heavy fragments(IMF's). There- 
fore, in case of FN's and LCP's transition energy tends 
to decrease at very stiff density dependence of symmetry 
energy, which is not in case of IMF's. The elliptical flow 
which is generated due to the compression in the partici- 
pant zone, tends to give different behaviour for the heavy 
fragments (generated from spectator zone) as compared to 
the light mass fragments (generated from the participant 
zone). 

In fig. 5, we show elliptical flow, (^2) at mid rapid- 
ity (|y| = < 0.1) for Z < 2 (left panel) and for 

protons(right panel) as a function of the incident energy. 
The rapidity cut is in accordance with the experimental 
findings. We compare the theoretical results with constant 
form of the symmetry energy and the density dependent 
symmetry energy with 7 = 0.66, with the experimental 
data extracted by INDRA, FOPI, and PLASTIC BALL 
Collaborations [4l[5j|6]. The elliptical flow shows a transi- 
tion from in-plane to out-of-plane with increase in incident 
energy. The elliptical flow is found to become more nega- 
tive with the increase in incident energy. This is because 
the mean field, which contributes to the formation of a 
rotating compound system, becomes less important, and 
the collective expansion process based on nucleon-nucleon 
scattering starts predominant. The out of plane emission 
decreases again toward the higher incident energies. This 
happens due to the faster movement of the spectator mat- 
ter after (W2) reaches a maximal negative value. The den- 
sity dependent symmetry energy tends to give the more 
squeeze out of the nuclear matter as compared to the con- 
stant form. 

For 7 = 0.66, the symmetry energy tends to play 
larger role across the high density region. The repulsion 
produced due to the symmetry energy tends to give the 



larger squeeze out flow, which is evolved due to the pres- 
sure generated in the participant zone. The trends ob- 
tained through the density dependent symmetry energy 
for 7 = 0.66 are in more agreement with the experimental 
data as compared to the constant form of the symmetry 
energy. However, no such effect is observed in case of the 
low incident energies. At low incident energies the density 
of the system is not so large to see a huge variation in the 
symmetry energy. As in the case of higher incident ener- 
gies of the order of 800-1000 MeV/nucleon, the density of 
the system is 2-3 times the normal nuclear matter density. 
As the mid rapidity region corresponds to the participant 
zone. The high density achieved in the participant zone, 
due to target and projectile collisions, tends to produce the 
effect of density dependent symmetry energy. The larger 
squeeze-out is also observed in case of protons. 



4 Conclusion 

We have checked the sensitivity of elliptical flow for the 
various forms of the symmetry energies. We found that 
the elliptical flow is highly sensitive towards the different 
forms of the density dependent symmetry energy. Transi- 
tion energy tends to give different behaviour for the lighter 
particles as compared to the heavy fragments. However, 
larger squeeze out is observed with the stiff density de- 
pendence of symmetry energy(7 = 0.66) as compared to 
the full symmetry energy strength i.e. 32 MeV. 
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